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The  hysteresis  motor  is  a  well-known  synchronous  motor  that  is  used  in 
special  small  power,  high  speed  applications.  Dynamic  modeling  and 
analysis  of  this  motor  is  more  complicated  than  permanent  magnet 
synchronous  motors  (PMSMs)  or  induction  motors  (IMs)  due  to  nonlinear 
behavior  of  rotor  magnetic  material.  Short  over-excitation  is  a  unique 
phenomenon  that  only  occurs  in  hysteresis  motor  in  which  the  terminal 
voltage  increase  at  synchronous  speed  for  a  short  duration,  and  then 
continuously  is  decrease  to  initial  value.  Therefore,  the  input  current  is 
reduced,  this  leads  to  more  power  factor  and  efficiency  enhancement.  Till 
now,  there  isn’t  any  analytic  dynamic  model  of  this  phenomenon.  In  this 
paper,  based  on  a  novel  dynamic  model  of  hysteresis  motor,  the  over¬ 
excitation  phenomenon  is  investigated  and  transient  performance  of  the 
motor  during  over-excitation  is  simulated  via  Simulink. 
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1.  INTRODUCTION 

The  hysteresis  motor  can  be  regarded  as  a  self-starting  synchronous  machine,  which  has  the  same 
structure  as  squirrel  cage  asynchronous  motor  with  exception  that  its  rotor  consist  of  a  ring  of  magnetic 
material  which  covers  a  non-magnetic  core.  It  is  often  used  for  small  power  application  that  needs  a  very 
smooth  torque  at  high  speed  such  as  gas  centrifuge  and  gyroscope  [1].  Relatively  low  starting  current,  simple 
and  strong  structure,  uniform  torque-speed  characteristic  and  also  low  power  factor,  low  efficiency  and 
hunting  phenomenon  are  some  of  hysteresis  motor  advantages  and  disadvantages  [2].  There  are  two  main 
structures  for  hysteresis  motors;  cylindrical  and  disk  types  as  shown  in  Figure  1.  Cylindrical  hysteresis 
motors  are  classified  in  two  types  of  circumferential-flux  and  radial-flux.  In  circumferential-flux  hysteresis 
motor,  the  rotor  ring  is  mounted  on  non-magnetic  material  as  support.  The  magnetic  field  lines  in  the  rotor 
are  mostly  circumferential  to  the  ring.  This  type  of  hysteresis  motor  is  mostly  used  in  industrial 
applications  [3]. 

Since  the  power  factor  of  hysteresis  motor  is  much  lower  than  that  of  an  induction  motor,  it  takes 
more  RMS  current  and  hence  large  capacity  condensers  are  installed  to  enhance  the  power  factor.  However, 
this  is  not  favorable  because  resonance  is  caused  by  the  presence  of  hysteresis  motor  inductance  and 
condenser  capacitance  and  efficiency  is  not  improved.  Rotors  showed  that  the  over -excitation  for  a  short 
period  of  the  hysteresis  motor  running  at  synchronous  speed  causes  the  reduction  of  the  stator  current  and  the 
increase  of  the  pullout  torque  at  the  same  time  [4] .  This  phenomenon  is  very  important  from  the  view  point  of 
high-efficiency  operation  and  has  been  discussed  by  some  researchers  [5]-[8].  They  have  investigated  the 
over-excitation,  but  couldn’t  model  this  phenomenon.  Lately,  over -excitation  has  been  experienced  via 
variable  frequency  inverters  and  various  voltage  patterns  have  been  applied  to  hysteresis  motor  to  gain  more 
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efficient  performance  for  hysteresis  motor  [9]-[l  1].  But  it  has  not  presented  any  dynamic  model  to  predict  the 
transient  performance  of  hysteresis  motor  during  over -excitation.  This  paper  investigates  the  transient 
behavior  of  the  motor  during  short  over-excitation  based  on  a  novel  dynamic  model  for  hysteresis  motor.  The 
paper  is  organized  as  follows:  In  section  2,  the  over-excitation  phenomenon  is  described  briefly.  Section  3 
proposes  a  novel  dynamic  model  of  hysteresis  motor.  In  section  4,  some  simulation  results  are  given  and  then 
conclusion  is  explored  in  section  5. 


Figure  1.  Structure  of  various  types  of  hysteresis  motor:  (a)  Circumferential -flux  (b)  Axial-flux 


2.  SHORT-DURATION  OVER-EXCITATION  PHENOMENON 

Short-duration  or  short  Over-excitation  is  an  effective  process  for  increasing  output  power, 
improving  power  factor  and  efficiency  and  reducing  input  current  [12].  Over -excitation  means  that  the  input 
voltage  Vi  applied  to  a  motor  operating  at  synchronous  speed  is  increased  continuously  up  to  nVi  (n  >  1)  and 
then  continuously  decreased  up  to  Vi.  Factor  n  is  called  over-excitation  factor. 

2.1.  Terminal  Voltage  Pattern  with  Short-Duration  Over-excitation 

Figure  2  shows  the  method  for  the  over-excitation  considered  here.  The  motor  is  accelerated  to  its 
synchronous  speed  at  a  stator  voltage  Vi.  After  synchronization,  the  stator  voltage  is  raised  to  a  higher  value 
nVi  for  a  short  period,  and  then  reduced  to  the  original  value  Vi.  It  is  assumed  that  the  load  torque  remains 
constant  during  this  change.  In  this  paper  a  circumferential -flux  type  hysteresis  motor  is  analyzed,  although 
the  following  method  can  also  be  extended  to  radial-flux  type  hysteresis  motor. 


Figure  2.  Terminal  voltage  pattern  [11] 


2.2.  Magnetic  Behavior  of  Rotor  Hysteresis  Material  before  and  after  Over -excitation 

Let  us  first  investigate  the  magnetic  state  of  rotor  before  over-excitation.  When  the  rotor  turns  at 
sub-synchronous  speed,  any  point  in  rotor  hysteresis  material  experiences  a  sinusoidal  variation  of  flux 


IJPEDS  Vol.  8,  No.  2,  June  2017  :  623  -  638 


IJPEDS 


ISSN:  2088-8694 


n  625 


density  B  in  time,  and  consequently  a  magnetic  hysteresis  loop  Lq  shown  in  Figure  3(a).  Other  points  in  the 
rotor  hysteresis  material  experience  exactly  the  same  hysteresis  loop,  displaced  in  time  phase.  Thus  the  loop 
not  only  gives  the  time  variation  of  magnetic  state  at  a  point  in  the  rotor,  but  also  gives  space  variation  of 
magnetic  state  around  rotor  at  a  particular  instant  of  time  [7,13].  The  area  of  this  loop  is  proportional  to 
output  torque  Tem-  When  the  motor  reaches  to  synchronous  speed,  magnetic  state  at  any  point  momentarily 
stops  changing.  For  example,  the  magnetic  state  at  a  point  P  in  Figure  3(b)  in  the  rotor  stops  at  the  location  of 
Po  on  the  loop  shown  in  Figure  3(a),  and  the  magnetic  state  at  another  point  stops  at  Qo  on  the  same  loop.  If 
the  load  torque  is  equal  to  output  torque  Tem,  motor  speed  will  not  change  further.  But,  if  the  load  torque  is 
less  than  Tgm,  the  motor  speed  will  increase  further  and  a  displacement  of  the  rotor  forward  with  respect  to 
the  rotating  flux  density  wave  will  occur.  So,  points  of  Pq-Qo  come  to  PrQi  that  means  the  operational  B-H 
loop  goes  to  other  loop  Fi.  Hence,  a  new  B-H  loop  is  created  that  its  area  is  smaller  than  loop  Fq.  It  means 
that  with  changing  load  torque,  the  B-H  loop  will  not  be  fixed.  The  same  behavior  can  be  explained  when  the 
load  torque  increases,  so  that  loop  Fi  will  be  wider  than  loop  Fq.  The  shape  of  Fi  is  different  from  the  major 
hysteresis  loop  created  by  rated  input  voltage  or  from  those  of  the  minor  hysteresis  loops  created  by  lower 
voltages.  Since  the  area  of  the  B-H  loop  is  proportional  to  the  output  torque  [1],  the  area  of  loop  Fi  is 
proportional  to  output  torque  witch  is  in  equilibrium  with  load  torque. 

Now,  let  us  investigate  the  change  in  the  magnetic  state  of  the  rotor  when  it  is  overexcited.  When 
the  stator  voltage  is  increased  from  Vi  to  nVi,  the  flux  density  in  the  rotor  also  increases.  Therefore,  the  points 
Pi,  Qi  on  loop  Fi  will  ascend  along  their  respective  minor  loops  toward  the  ascending  branch  of  the  major 
loop  Fq  in  Figure  4(a).  The  points  that  meet  the  ascending  branch  of  loop  Fq  will  go  up  further  along  it  the 
points  that  pass  the  tip  of  the  major  loop  Fq  will  ascend  along  the  normal  magnetization  curve  of  the  rotor 
hysteresis  material.  For  example,  the  points  PrQi  will  come  to  P2-Q2,  respectively.  Thus,  the  B-H  relation 
around  the  rotor  is  given  by  a  closed  loop  F2  in  Figure  4(a).  It  should  be  noted  that  load  torque  remains 
unchanged  during  mentioned  process.  This  means  that  the  area  of  loop  F2  must  be  equal  to  that  of  loop  Fi.  In 
other  words,  amplitude  and  phase  of  flux  density  distribution  in  the  rotor  are  such  that  the  areas  of  both  loops 
are  equal.  Such  flux  density  distribution  and  therefore  loop  F2,  can  be  obtained  from  an  iterative  calculation. 
When  stator  voltage  is  reduced  from  nVi  to  original  value  Vi,  the  flux  density  in  rotor  now  decreases. 
Therefore,  for  example,  points  P2  and  Q2  on  loop  F2  descend  along  the  descending  branch  of  the  respective 
major  loops  whose  tips  are  P2  and  Q2,  and  come  to  P3  and  Q3,  respectively,  as  shown  in  Figure  4(a).  Thus,  the 
B-H  relation  around  the  rotor  is  given  by  a  closed  loop  F3.  The  area  of  loop  F3  must  be  equal  to  that  of  loop 
F2,  because  the  load  torque  is  constant.  If  the  load  torque  is  increased,  there  will  be  a  decrease  of  speed  and  a 
displacement  of  the  rotor  backward  with  respect  to  rotating  flux.  Then,  for  example,  the  flux  density  at  point 
P  decreases  to  P4  and  the  flux  density  at  point  Q  increase  to  Q4  as  shown  in  Figure  4(b).  As  a  result,  the  B-H 
relation  given  by  loop  F4  in  Figure  4(a)  can  be  obtained,  the  area  of  which  is  corresponding  to  the  pull-out 
torque  after  over-excitation  [7]. 


Figure  3.  B-H  loops  and  flux  density  waves  for  two  Figure  4.  B-H  loops  and  flux  density  waves  during 
different  loads  [7,13],  (FI  loop  corresponds  to  the  over-excitation  process  [7] 

lower  load) 
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3.  DYNAMIC  MODEL  of  HYSTERESIS  MOTOR 

The  common  point  in  all  mentioned  steady-state  equivalent  circuit  models  is  that  the  rotor's 
hysteresis  loop  is  considered  as  hysteresis  resistance  Rh  in  series  with  reactance  Xh  in  sub  synchronous  mode 
as  shown  in  Figure  5.  In  this  section  a  new  dynamic  model  for  hysteresis  motor  with  incorporating  the 
change  of  load  torque  and  stator  voltage  to  identify  the  operational  B-H  loop  of  rotor  is  developed  [14]. 


Figure  5.  Equivalent  circuit  of  hysteresis  motor  at 
sub  synchronous  mode  [15] 


3.1.  Calculation  of  Rotor’s  Equivalent  Circuit  Parameters 

Hysteresis  motor  used  in  this  paper  is  a  super-high  speed  circumferential -flux  type  motor,  3 -phases, 
star-connected,  with  given  parameters  in  Table  1.  To  develop  dynamic  model  of  hysteresis  motor,  the  value 
of  rotor’s  parameters  due  to  operational  hysteresis  motor  must  be  calculated.  For  this  purpose  B-H  curves  of 
rotor  material  are  firstly  measured  by  some  standard  experiments  [16].  Employed  magnetic  material  used  in 
this  study  is  Vicalloy  2  with  typical  B-H  loops  shown  in  Figure  6.  Maximum  achieved  output  hysteresis 
torque  in  steady  state  is  corresponding  to  maximum  applied  voltage  to  stator  that  is  proportional  to  major 
hysteresis  loop.  For  simplicity  of  calculation,  the  fundamental  harmonic  of  magnetic  field  intensity  {H) 
and  field  density  {B)  are  just  considered  that  means  the  hysteresis  loops  are  considered  as  elliptical. 


Table  1.  Rated  specifications  and  parameters  of  used  circumferential-flux  hysteresis  motor 


Symbol 

Quantity 

Value 

Dimension 

^rated 

line  rated  voltage 

230 

V 

T 

^  rated 

rated  torque 

0.01 

N.m 

f rated 

rated  frequency 

1000 

Hz 

CO 

supply  angular  frequency 

271x1000 

rad/sec 

m 

number  of  phase 

3 

P 

number  of  poles 

2 

J 

shaft  inertia  moment 

3x10^“* 

kg.m^ 

Rs 

stator  resistance 

16.4 

Q/ph 

Rc 

stator  core  loss  equivalent  resistance 

10580 

Q/ph 

Rh 

rated  rotor  hysteresis  resistance 

300 

Q/ph 

Re 

eddy  current  resistance  of  rotor 

223 

Q/ph 

X,s 

stator  leakage  reactance 

78 

Q/ph 

rated  magnetizing  reactance 

400 

Q/ph 

Xh 

rated  rotor  hysteresis  reactance 

170 

Q/ph 

h 

axial  length  of  hysteresis  ring 

25 

mm 

g 

air  gap  length 

1 

mm 

mean  radius  of  air  gap 

24.45 

mm 

A 

hysteresis  lag  angle 

60^ 

p 

electrical  resistivity  of  hysteresis  ring 

650 

pQ-mm 

Pr 

relative  permeability  of  hysteresis  ring 

20 

Po 

permeability  of  free  space 

471X10'^ 

H/m 
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Figure  6.  B-H  characteristics  for  a  typical  Vicalloy  magnetic 
material  used  in  hysteresis  motor  [17] 


The  flowchart  for  calculation  of  the  motor’s  parameters  in  air  gap  and  rotor  circuit  is  shown  in 
Figure  7.  Firstly,  the  operational  B-H  loop  or  corresponding  lag  angle  ^  due  to  applied  stator  voltage  is 
identified  and  then  is  modified  based  on  applied  load  torque  using  equivalent  circuit  shown  in  Figure  5. 

In  this  circuit,  K  +JX,  is  stator  leakage  impedance,  II  jX  ^  is  air  gap  impedance,  and 
(R,+jX,)\\R^  / s  is  equivalent  rotor  impedance.  The  air  gap  reactance  and  rotor  impedance  are  expressed 
as  [18]: 


X 


m 


KP^g 


(1) 


R, 


_mco(K^N  , 


siny^ 


(2) 


_mco[K^X  ^ 


Tt^rl 


COSy0 


(3) 


To  calculate  the  eddy  current  resistance  in  hysteresis  material  of  the  rotor,  it  should  be  noted  the  flux 
flow  in  rotor  is  circumferential  and  the  eddy  current  vector  is  perpendicular  to  the  flux.  So,  eddy  current 
resistance  can  be  obtained  from: 


The  hysteresis  ring  volume  Vy  and  cross  sectional  area  Ah  in  (2)-(4)  are  calculated  from: 


(4) 

(5) 

(6) 
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Figure  7.  Flowchart  for  calculation  of  the  rotor's  equivalent  circuit  parameters  [18] 


To  identify  the  hysteresis  lag  angle  ^  due  to  operational  stator  voltage  firstly  an  initial  B-H 
curve  (or  initial  and  5^)  is  chosen.  Suppose  the  fundamental  harmonic  of  magnetic  field  intensity  in  the 
rotor  is  defined  as: 


H  =  H^co^cot  =  {B^I lu)co^cot  ^  ’ 

So,  flux  density  in  rotor  can  be  expressed  as: 

Corresponding  to  measured  B-H  loops,  parameters  of  Cl->b  can  be  calculated  for  any  given  and  created 
Byn.  For  an  elliptical  B-H  loop  with  area  of  Wh,  the  parameter  of  h  is  related  to  Wh  through: 


(9) 


and  coefficient  ‘a'  is  derived  from: 


(10) 


Therefore,  the  hysteresis  lag  angle  ^  is  obtained  by: 


y^  =  tan  \bla) 


(11) 
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which  it  varies  from  40°  to  70°  for  semi-hard  magnetic  material  used  in  this  study.  Now  as  shown  in  Figure  8, 
for  obtained  lag  angle  P,  the  equivalent  circuit’s  parameters  are  calculated.  To  check  the  validity  of  P, 
computed  input  voltage  of  stator  is  compared  with  real  stator  voltage.  To  compute  the  stator  voltage  Vy, 
induced  voltage  in  air  gap  is  calculated  from: 


Figure  8.  Equivalent  circuit  for  dynamic  model  of  hysteresis  motor  in  rotating  dq  reference  frame 


where  (j^g  represents  the  air  gap  flux  per  pole,  and  is  obtained  for  a  cylindrical  machine  with  selected  value 
of  by: 


<4 

The  stator  input  current  from  steady-state  equivalent  circuit  shown  in  Figure  5  can  be  obtained  from: 


L  =. 


^g  ^g  n 

—^  +  — —  +  —^cosS 

Rc 


+ 


^g  ^g  ■  n 

m  J 


^2 


(13) 


(14) 


that  the  stator  phase  voltage  is  considered  as  reference  phasor  and  Og  is  the  phase  difference  between  E 


and  Is  as: 


d„  =tan~ 


1  sin  6 

- ^ - — 


1  1  cos P 

\K  ^  J 


(15) 


and  the  impedance  due  to  hysteresis  material  is: 

Therefore,  the  input  voltage  is  calculated  by: 

V in  =  ^l{EgCos0g  +ZJsCos0^  f  +(^Eg  sin  0^  +2^1^  sin  0^  "f 


(16) 


(17) 


where  are  the  stator's  impedance  and  argument  as: 
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Z,=Jr/7x^, 

The  error  between  calculated  input  voltage  by  (17)  and  real  stator  voltage  should  be  less  than 
Si  =  0.01  y  .  Otherwise  the  chosen  B-H  curve  is  changed  by  modifying  of  by  a  proper  step  size,  and  the 
calculations  are  repeated  until  the  error  falls  into  Si  tolerance.  As  mentioned  earlier,  when  the  hysteresis 
motor  reaches  synchronous  speed,  developed  torque  Tem  is  equal  to  load  torque  Ti.  The  steady-state 
developed  torque  in  terms  of  hysteresis  ring’s  volume,  operational  B-H  loop  properties  and  lag  angle  P  can  be 
expressed  as  [19]: 


fflP 

Tern  =-^^(  Wh  P 


(19) 


It  means  that  if  under  constant  stator  voltage,  the  load  torque  varies,  operational  B-H  loop  changes 
without  changing  value  of  and  and  so  the  lag  angle  P  varies  as  shown  in  Figure  3.  Therefore,  in  final 
stage  of  computation,  developed  torque  through  (6)  is  compared  with  the  load  torque  and  lag  angle  P  is 
iterated  by  a  proper  step  size  until  the  torque  difference  falls  to  desired  tolerance  S2  =0.0005  N.m  as 
shown  in  Figure  7.  The  computed  lag  angle  P  is  then  used  to  update  the  rotor's  equivalent  circuit  parameters. 

Hysteresis  motors  use  the  semi-hard  magnetic  materials  that  unlike  to  induction  motors  induce 
considerable  back-EMF  voltage  in  stator.  For  dynamic  modeling  of  hysteresis  motor,  it  is  necessary  to 
consider  an  induced  voltage  source  in  rotor  circuit.  This  equivalent  induced  voltage  for  hysteresis  material 
can  be  obtained  from  Figure  5  as: 

Ep  —  g  cos  0g  +{X  -\-X ^  cos 


Where 


61,  =tcin 


XJsCosOg  +(X;,  +X^)lf^cosP 
XmlsSinOg  +X^)l^sinP 


(21) 


Calculated  Eh  is  a  phasor  corresponding  to  phase  ‘a’.  It  has  to  be  expressed  in  time  domain  for  three  phases 
and  then  transferred  to  E^h  and  Eqh  in  dq  reference  frame.  These  voltages  at  synchronous  condition  are 
constant. 

3.2.  Mathematical  Model  and  Equivalent  Circuit  of  Hysteresis  Motor  in  Rotating  dq  Reference  Frame 

The  hysteresis  motor  voltage  equations  in  synchronously  rotating  dq  reference  frame  can  be 
written  as: 


^dq  -^  '^dq  + 


dt 


+  E, 


dq 


where. 


1 

1 _ 

1 

_ 1 

_ 1 

II 

o 

II 

^qs 

^dr 

II 

hr 

II 

^qs 

^dh 

1 

o 

II 

^qr 

1 

1 _ 

^qh 

(22) 


(23) 
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and, 


R 


0  0  0 

0  /?,  0  0 

0  0  0 
0  0  0  /?,. 


(24) 


where  the  rotor  resistance  is  calculated  from: 

R.  =  ^ 


l/K+s/R, 


(25) 


and  ‘5*’  denotes  the  slip.  The  stator  and  rotor  fluxes  in  (22)  are: 


Ajq  - 


where  the  inductances  are  calculated  from: 
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The  magnetizing  and  rotor  inductances  can  be  calculated  from  (1)  and  (3).  Two  voltages 

in  (23)  are  defined  by: 


^■'ds  ^r^qs  ^ 


^qs  ~ 


r^ds 


(28) 


Also  Edh  and  Eqh  are  dq  components  of  induced  voltage  due  to  residual  magnetizing  action  of  hysteresis 
material  on  the  rotor  that  are  obtained  using  (20)  as  follows: 


-^dh 

-^qh 


■y40r)_ 


Ef^  sin  cot 

Efj  sin( ojt-HQR ) 

Ej^  sin(  QX  +  Vl(f  ) 

That  the  following  Park  matrix  is  used  for  abc  to  dq  transformation  [20]-[22]: 


(29) 


[W^r)]  =  | 


cosO^  cos( 0J.  -120°  )  cos( 6^  4-120°  ) 
sin  6^  sin(  9^  -120°  /  sin(  9^  4-120°  ) 
11  1 
.2  2  2  . 


(30) 


Developed  electromagnetic  torque  can  be  expressed  in  term  of  flux  and  current  as: 

3  P 


(31) 


=■ 


2  2 


(  ^ds^qs  ^qs^ds  ^ 


and  the  mechanical  speed  is  obtained  from: 
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dco^(  t ) 
dt 


=  T 

^em 


-Ti 


(32) 


The  foregoing  equations  can  be  summarized  in  the  equivalent  circuit  model  as  shown  in  Figure  8. 


4.  SIMULATION  RESULTS 

Some  simulations  are  explored  for  a  hysteresis  motor  with  parameters  summarized  in  Table  1. 
It  should  be  noted  that  employed  hysteresis  motor  in  this  study  is  a  super  high-speed  motor  with  rated  speed 
60,000  rpm  in  which  its  start-up  procedure  from  zero  to  final  speed  takes  4200  sec  with  a  slow  acceleration 
rate.  Simulation  of  this  duration  takes  very  much  computation  and  needs  a  lot  of  RAM  and  memory. 
Therefore,  to  shorten  the  simulation  time  to  2  seconds,  we  are  forced  to  decrease  the  inertia  moment 
adequately.  Before  simulation  of  over-excitation  phenomenon,  developed  dynamic  model  is  simulated  to 
confirm  the  theoretical  aspects  of  hysteresis  motor.  To  have  better  realization,  some  of  quantities  are  given  in 
per  unit  scale.  The  predominant  load  of  employed  hysteresis  motor  is  rotational  friction  as  -col  that  the 

value  of  Bf  is  1  in  per  unit  scale.  Although  an  extra  constant  load  at  synchronous  speed  is  applied  to 
hysteresis  motor. 

Figure  9  to  Figure  12  show  the  simulation  results  of  hysteresis  motor  during  start-up  with  open-loop 
V/f  scalar  control  scheme.  In  addition  to  rotational  friction  load,  a  constant  load  is  applied.  Figure  9  shows 
the  applied  reference  frequency  and  applied  stator  phase  voltage.  In  figure  10,  the  speed  and  torque  response 
of  motor  are  depicted  during  startup.  The  speed  reaches  synchronous  speed  under  friction  and  constant  load. 
The  speed  tracking  error  during  start-up  is  at  most  0.5%  and  becomes  zero  at  steady-state.  Moreover,  when  a 
constant  load  torque  0.5  pu  is  applied  at  time  1.5  sec,  the  motor  can  reach  to  synchronous  speed  after  some 
oscillations.  The  RMS  current,  input  power  and  power  factor  waveforms  are  shown  in  Figure  11.  The  current 
under  rated  torque  1  pu  at  steady-state  (t=l.2sec)  is  0.19  A  (close  to  actual  current  of  the  motor  in 
experiment).  When  the  load  torque  increases  to  1.5  pu,  the  current  increases  small  value  and  reaches  to  0.23 
A.  It  is  accordance  to  behavior  of  an  actual  hysteresis  motor.  The  input  power  of  the  motor  at  steady  state  is 
42  W  includes  of  output  power  and  stator's  copper  and  core  losses.  With  increasing  the  load  torque,  the 
power  also  increases  50%.  Also  the  power  factor  firstly  is  0.54  and  after  1.5  sec,  increases  to  0.72. 
In  Figure  12  the  variations  of  hysteresis  lag  angle  and  rotor  parameters  are  displayed.  The  angle  p  is 
increasing  as  well  as  the  stator  voltage  increases  during  V/f  and  after  1  sec  it  converges  to  55°.  Also,  at  time 
1.5  sec  and  after  increasing  the  load  torque,  it  converges  to  higher  value  60°.  It  means  that  during  startup  or 
load  changing,  the  operational  B-H  loop  continuously  changes  and  finds  corresponding  B-H  loop  after  1  sec 
when  the  voltage  and  speed  become  constant.  Also,  the  values  of  are  changing  with >0 

through  (8)-(9). 


Syn  chronou  s  frequency 


1000 

"nT 

X 

500 

0 

0„5  1  1 3 

Stator  phase  voltages 


200  - 
o  - 
-200 


OJZ  0.4  0.6  o.a 

time  (Sec) 


Figure  9.  Reference  frequency  and  applied  phase  voltage  to  hysteresis  motor  during  open-loop  V/f 
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Figure  10.  Speed  and  torque  response  of  hysteresis  motor  during  open-loop  V/f  scheme  start-up 
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Figure  11.  RMS  current,  input  power  and  power  factor  of  hysteresis  motor  during  open -loop  V/f 

scheme  start-up 
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Figure  12.  Variations  of  hysteresis  lag  angle  (3  and  Rj^  ,Xj^  of  rotor’s  circuit  model 
Now,  the  short-duration  over-excitation  is  simulated  based  on  developed  dynamic  model.  Figure  13 
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shows  the  current  and  power  factor  of  motor  for  over-excitation  factor  (n=1.25).  Input  current  decreases  after 
reduction  input  voltage  from  0.19  A  to  0.13  A.  It  means  the  cupper  loss  of  stator  reduces  4  times  than  before 
over-excitation.  On  the  other  hand,  power  factor  increases  from  0.54  to  0.77  that  means  a  50%  enhancement. 


RMS  of  phase  voltage 


X:0.5 


Power  factor 


Figure  13.  Input  voltage,  current,  power  factor  before  and  after  over-excitation  with  n=1.25 


Figure  14  shows  the  variations  of  lag  angle  (3  and  rotor’s  impedance  components  due  to  hysteresis 
material  of  rotor.  During  over-excitation,  lag  angle  decreases  and  then  it  increases.  The  variations  of  Rh 
and  Xh  obey  (2)-(3). 


Hysteresis  lag  angle 
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Figure  14.  Lag  angle  (3,  Rh  and  Xh  before  and  after  over-excitation  with  n=1.25 


Figure  15  shows  the  current  and  power  factor  of  motor  for  n=l.l.  The  RMS  value  of  current 
decreases  form  0.19  A  to  0.16  A,  and  the  power  factor  increases  from  0.54  to  0.64.  Comparing  of  the  results 
shown  in  Figure  13  and  Figure  15  proves  that  the  larger  increase  of  terminal  voltage  pattern  leads  to  better 
performance  and  higher  power  factor  and  efficiency.  Using  developed  dynamic  model  of  hysteresis  motor,  it 
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is  possible  to  predict  the  best  choice  for  order  “n”  for  over-excitation  process. 
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Figure  15.  Input  voltage,  current,  power  factor  before  and  after  over-excitation  with  n=l.l 


Figure  16  shows  the  result  of  over-excitation  process  with  step  change  of  terminal  voltage 
and  n=1.25.  The  enhancement  of  power  factor  and  reduction  in  RMS  current  is  identical  to  ramp  variation  in 
terminal  voltage. 
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Figure  16.  Input  voltage,  current,  power  factor  before  and  after  over-excitation  for  step  change  of 

voltage  with  n=1.25 
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However  as  shown  in  Figure  17,  the  fluctuation  of  speed  for  step  change  of  the  voltage  is  higher 
than  ramp  change. 
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Figure  17.  Changes  of  slip  and  voltage  before  and  after  over-excitation  for  (a)  ramp  change  (b)  step  change 

of  terminal  voltage  with  n=1.25 


5.  CONCLUSION 

Short  duration  over-excitation  increases  the  efficiency,  power  factor  and  pull-out  torque  of 
hysteresis  motor  and  makes  it  possible  to  operate  in  efficient  condition.  Determination  of  voltage  amplitude 
and  time  duration  in  over-excitation  process  is  essential  to  attain  the  best  operational  point.  In  this  paper,  the 
short  over-excitation  phenomenon  of  a  hysteresis  motor  has  been  modeled  and  simulated  for  the  first  time. 
On  this  way,  a  novel  dynamic  model  of  hysteresis  motor  in  rotation  dq  reference  frame  was  presented.  In 
developed  model,  the  operational  B-H  curve  of  the  motor  is  considered  during  change  of  input  voltage  or 
load  torque.  Simulation  results  show  the  effectiveness  of  developed  model  to  predict  the  motor  performance 
during  over-excitation.  The  results  confirm  that  with  higher  change  of  ststor  voltage,  the  stator  current 
reduces  more  and  so  the  power  factor  increases.  Moreover,  step  change  in  input  voltage  leasd  to  more 
hunting  of  speed.  By  using  developed  model  it  is  possible  to  evaluate  the  effect  of  motor  parameters,  B-H 
characteristics  of  rotor  material,  magnitude  of  input  voltage,  time  duration  and  acceleration/deceleration  of 
voltage  variations  during  short-over-excitation  process. 
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